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Abstract 

A technique is presented for treating strongly nonstationary and transient processes 
in optics, permitting one to take into account both types of competing with each other 
effects, quantum as well as coherent. The main equations for describing the interplay 
between these two kinds of effects are derived. The possibility of influencing coherent 
optical phenomena by preparing special quantum states of matter is discussed. 
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1. Introduction 

Coherent and quantum phenomena are often treated as antagonists, since the former are, 
to some extent, close to classical ones 1 . When describing one of these phenomena, one usually 
employs rather different approximations. For example, a very common way of considering 
coherent phenomena is by invoking the semiclassical approximation, which makes it possible 
to give a transparent picture of such phenomena. But the semiclassical approximation kills 
all quantum correlations, eliminating by this all quantum effects. Such a neglection of the 
latter may be admissible if the studied coherent process is stationary or a strong coherence 
is imposed on the system by an external field. However, if one investigates a self-organized 
growth of coherence and when the coherent phenomena are transient or intermittent, then 
quantum effects may essentially influence the features of coherent phenomena. In that case, it 
is necessary to take account of both types of effects. For this purpose, one may resort to the 
consideration of higher-order correlation functions. But then, the evolution equations become 
so much complicated that essentially nonstationary states are hardly treatable. 

In this report, a method is presented, which combines the transparency of the semiclassical 
approximation with the possibility of taking account of quantum effects. The idea of the 
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method reminds, to some extent, quantization of soliton solutions in quantum field theory . 
But the most important concept of the method is the notice that quantum and coherent effects 
occur on rather different scales. Quantum effects are usually connected with short-range fast 
fluctuations, while coherent phenomena are mainly long-range in space and slow in time, as 
compared to the quantum ones. The occurrence of different spatio-temporal scales allows the 
development of Scale Separation Approach . The general idea of separating different scales 
is, of course, known. However, its realization, as applied to optics in the report below is new. 
The pivotal novel technique is Quantization of Local Fluctuations. 

2. Resonant Atoms in Matter 

The system of resonant atoms inside matter is described by the Hamiltonian 

H = H a + Hf + H m + H a f + H m f . (f ) 

Here, the first term 

tf a = X> (5 + -J (2) 

is the Hamiltonian of resonant atoms, with transition frequency uq, where Sf are pseudospin 
operators. The field Hamiltonian is 

^=^/( E2 + B2 )*' P) 

The Gaussian system of units is used, with setting H = 1. Electric and magnetic fields are 
expressed through the vector potential, 

E = --^, B = VxA. 

c at 

The Coulomb gauge V • A = will be employed. This calibration is convenient for eliminating 
field variables in the following evolution operations. The Hamiltonian H m models the material 
incorporating the atoms. The atom-field interaction is presented by the dipole Hamiltonian 

N /I \ 
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in which Aj = A(i"j, i) is the vector-potential operator, E j = E (rj,t) is a classical external 
field. The transition current and transition dipole operators are 

J, = iu Q (dSf - d*S~) , P, = dS+ + d*S~ , (5) 

where d is a transition dipole and Sf = Sf ± iSf. Finally, the matter-field interaction is 
described by the Hamiltonian 

H m f = ~ -Jjm(r,t) ■ A(r,t) dr , (6) 
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in which j m (r, t) is a local density of current in matter. 

The Heisenberg equations of motion for the pseudospin operators yield the equation 



dt 



= -icuoSr + 2 (A; d • Aj — id • E w ) S* 



(7) 



for the lowering operator and the Hermitian conjugated equation for the rising operator, with 
k = uq/c. The equation for the population-difference operator is 



dS[ 
dt 



(k d • Aj - id • E 0l ) 5+ - (fc d* • A, + *d* • E 0l ) S~ . 



(8) 



In deriving an equation for the vector potential, one has to use the commutation relations 

~E a (r,t), Af(r',t)] = AmcS afi (r - r') , 



where 



<W(r) 



5, 



a(3 



k a k^ 



e ik r dk 



is the so-called transverse delta-function 1 . These commutation relations explicitly take into 
account the Coulomb gauge condition. The transverse delta-function can be presented in other 
forms, as 

or as 

3 6{T) ' 



<W(r) 



where n a = r a /r. As a result, for the vector potential one has 



1 <9 2 \ 
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^-?S»J A = -T J - 

with the density of current 

E Mr " r D^f(*) + / - r')j£(r', t ) rfr' 
.i=i J 

The solution to Eq. (9) reads 



r(r,o = E 

/3 Li=l 



(9) 



(10) 



A(r,t) = A mc (r,t) + - /j r',t 



(11) 



with A vac being the solution of the uniform equation related to Eq. (9). Note that the dipolar 
part of the transverse delta-function, being averaged over space, yileds zero, that is 



This anisotropic dipolar part, for simplicity, can be omitted. Then the vector potential (11) 
takes the form 

A = A vac + A. ra( i + A. ma t , (12) 

in which the term 

N 2 ( r ij 



is due to radiating atoms and the potential 

A mat (r,t) = ^-J j m (r',t 



2 r . ( , 1 „\ dr' 

— r — i 



r'l 



is caused by matter currents. Here = | r^- j , = — rj. The summation over j = 1,2, ... ,N 
does not include the term with j = i corresponding to self-action, which will be taken into 
account by incorporating the level and line widths in the evolution equations. 

It is important to notice that all processes are defined for t > 0, while in the expressions 
above there appears the dependence of operators on the difference t — t' , which may be nega- 
tive. Therefore it is necessary to complete the definition of quasispin operators by adding the 
retardation condition 

Sr(t)=0 (f<0). (13) 

The dependence of S^(t — t f ) on the retarded time can also be simplified by taking into account 
that the interaction of radiation with atoms is to be essentially weaker than intra-atomic inter- 
actions. In the other case, the very notion of resonant atoms, having a well-defined transition 
frequency, would not have sense. Then, Eq. (7) shows that S~(t) ~ exp(— iu t). Combining 
this with the retardation condition (13) gives 

S~(t - t') = Q(t - t')S~(t) exp^ot') . (14) 

Hence for the typical dependence on the retarded variables, one has 

(t - ^) = © (ct - r^Sj (t) exp (ikonj ) . (15) 

Such a simplification, combining the Born approximation with the retardation condition (13), 
can be called the Retarded Born Approximation. 

It is worth mentioning that it is sufficient to understand all operator equations in the weak 
sense, as the equations for appropriate matrix elements or averages. 

3. Separation of Local Fields 

The influence of vacuum and matter on an atom enters the evolution equations through the 
expression 

£(r,i) = 2k d-(A vac + A mat ) . (16) 
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Actually, electromagnetic vacuum and matter form an effective vacuum acting on an atom by 
means of the local field (16). Thus, one may say that there are two types of operator variables, 
the quasispin variables and the local field (16). An operator F, being a function of these 
variables, can, for brevity, be written as F(S,£), with S denoting the quasispin set {Sf} and 
£, the set {£(r, t)} of local fields at different points r. According to the existence of two types 
of variables, one may define two kinds of statistical averages over the related variables, either 
quasispin or local-field ones. Such partial averages are defined by means of the corresponding 
restricted traces, either 

<F> = Tr s pF(S,0, (17) 
with p being a statistical operator, or 

«F» = Tr£pF(S,0- (18) 

When averaging over quasispin variables, one may assume the validity of the mean-field 
decoupling 

<S?Sf> = <S?xSf> (i^j), (19) 

since atoms interact with each other through effective long-range forces. This decoupling is not 
equivalent to the semiclassical approximation as the local-field variables have not been involved. 
Then the atomic average < S'f > is, actually, an operator function of local fields. The latter 
describe local quantum fluctuations. 

For what follows, it is convenient to use the notation S"(t) = 5 ,Q? (r i , t). Also, the summation 
over atoms can be replaced by the integration over the sample, according to the rule 

g^'/ dr i" s v)- (20) 

The behaviour of atoms can be described by the following averages. For an atom at the 
point r, one may write a transition function 

u(r,t) = 2 < S~(r,t) > , (21) 
w(r,t) = u + (r,t)u(r,t) , (22) 



intensity of coherence 
and population difference 



s(r,t) = 2 < S z (r,t) > . (23) 
The effective force acting on an atom is given by the sum 

f(r,t) = f (T,t) + f rad (T,t)+Z(r,t), (24) 

in which 

/o(r,*) = -2id.E (r,*) (25) 
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is due to a classical external field, 

f rad (r, t)=2k < d • A rad (r, t) > (26) 

is caused by the radiation of other atoms, and £(r, t) is the local field (16). The radiation force 
(26) explicitly writes 

f rad (r, t) = -i lQ p J [G(r - r', t) u(r', t) - e 2 d G*(r - r', t) u+(v', t)] dr' , (27) 
where the transfer function is 

G(r,t) = e(ct-r)^— 

and the notation 

2 

7o = g ^o^o > r = l r l > d = d o e d , d = |d| 

is used. The quantity 70 is a natural half-width. 

In this way, for the functions u = u(r, t), w = w(r, t), and s = s(r, t), we derive the evolution 
equations 

— = -(iuj + j 2 )u + fs , (28) 



dt 

= -2 l2 w + (u + f + f+u) s , (29) 

! = -^(W + / + «)-7i(*-C), (30) 

in which 71 is a longitudinal relaxation parameter, 72 is a transverse attenuation parameter, 
and ( G [—1, 1] is a stationary population difference for a single atom. 

4. Sample of Cylindrical Shape 

Let the sample have typical for lasers cylindrical shape. The axis of the cylinder is along 
the z-axis, which is distinguished by the propagating field 

E (r, t) = - E 1 e i(fc *- fcrt) + - E*e^ ( ^"^ , (31) 

where the frequency u = kc is close to the transition frequency u> , so that the resonance 
condition 

— <&1 (A = ou - u ) (32) 
^0 
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holds. The wavelength A = 2-kc/uj is small as compared to the radius, R, and length, L, of the 
cylinder, 

A « 1 , A « i . (33) 

Assuming the absence of sharp transverse structures, one may employ the single-mode approx- 
imation 

u(r, t) = u(t)e ikz , w(r, t) = w(t) , s(r, i) = s(f) . (34) 

The absence of sharp transverse nonuniformity implies that there exists the main propagating 
mode selected by the field (31). Expressions (34) are to be substituted in Eqs. (28) to (30). 
Equation (28) is multiplied by e~ lkz and then all equations are averaged over space. The 
following notation will be used for an effective force 

/!(*) = -id- Exe"*-* + (35) 

in which 

Z(t) = ±Jar,t)e-* kz dT. (36) 
Let us also introduce the coupling functions 

a(t) S 7o p jeict- r) Bin( ^; kz) dr , «<) S w / B(ct - r) ""^ fc » * . (37) 

Finally, Eqs. (28) to (30) are transformed into the ordinary differential equations 

du 

— = - [i(u + (3s) + 72 - as] u + /is , (38) 
^ = _2( 72 - as)w + («+/! + f+u) s , (39) 

§ = -<*» ~ \ {u + fi + A + u) - 7i(« - ■ (40) 

Although, one should remember that these are, actually, operator equations with respect to 
the quantum variable £(t). 

5. Method of Stochastic Averaging 

To simplify further the evolution equations (38) to (40), one can take into consideration the 
existence of several small parameters, such as 

— <1, — <1, — <1. (41) 

LOq LOq UJq 

The amplitude of the external field (31) is assumed to be small, so that 

1^ < 1 , v x = d ■ Ei . (42) 
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Also, the local quantum field £ is treated as weak, in the sense that its first and second moments 
are proportional to values much smaller than c^o- Then the multiscale averaging technique can 
be generalized ~ to stochastic and operator equations, as Eqs. (38) to (40). The occurrence of 
the above small parameters shows that the functions w(t) and s(t) are temporal quasi-invariants 
with respect to the fast function u(t). 

Equation (38) for the fast function, with w and s being quasi-invariants, can be solved. 
To this end, let us introduce the collective width and collective frequency by the corresponding 
expressions 

r = 7 2 - as , n = cu + f3s , (43) 
and also, let us define the dynamical detuning 

5 = uj-n = A-(3s . (44) 

The solution of Eq. (38) reads 

( Uo - jv) e -(*+n* + g-*-* + 8 f ^ )e -an + m-t') dt , (45) 

V o + iT J o + iT Jo 



u 



To simplify the following formulas, it is convenient to choose the phase of Ex so that to eliminate 
the dependence on the transverse initial value u . For this purpose, the phase of E 1 is taken so 
that u^d ■ Ex be real, which writes 

u%v x = u vl . (46) 

The solution (45) is to be substituted into Eqs. (39) and (40), whose right-hand sides are 
to be averaged according to the prescription 



<F> = lim i / T <F(£,t)> dt, 

T-^oo 1 Jo 



where the integration over time does not touch quasi-invariants. The quantum field £(£) is 
centered so that 

« m » = . (47) 

To make the resulting expressions less cumbersome, let us consider the case of small detuning 
\5\ < |T|, when 5 can be omitted in the phase dependence, though should be kept in de- 
nominators to avoid spurious poles. In the process of the averaging, one obtains the effective 
attenuation 

f=j^(i-e- r, ) +r 3. m 

where the first term is due to the classical external field (31) and the quantum attenuation 

T 3 = Re lim i F dt f < t{*W) > e^ in+r ^-^ dt' (49) 

T^oo T Jo Jo 

appears because of the action of the quantum field £(£). 
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Equations (39) and (40) reduce to the evolution equations 

^ = -2( 72 - as)w + 2Ts 2 , ^ = -aw - Ts - 7l (s - C) , (50) 
at at 

describing the coherent guiding centers. In order to make these equations complete, it is 
necessary to define the quantum correlation function <C £ + (£)£(£') 3> entering the quantum 
attenuation (49). 

6. Examples of Quantum Attenuation 

To make it clear how the quantum attenuation (49) can be calculated, let us give some 
examples of defining the quantum variable £(£). 

The first simple case could be by considering the variable £(£) as random, associated to 
infrared noise characterized by 

«£ + W£(0» = 73 2 - (51) 

Then Eq. (49) yields 

r « = TFTF • (52) 

The opposite case would be to treat as a stochastic variable representing white noise, 
with 

« e(t)m » = 2r 3 s(t - o , (53) 

which results in the identity = IV 

A more elaborate modelling of the effective quantum fluctuations is by a system of oscilla- 
tors, yielding 

eW = E7K)(v^ + ^) , (54) 

i 

where u q = u^ q > 0. For the Bose operators b q and b q , statistical averaging gives 

< b q > = , < b\b p > = n q 5 qp , < 6 9 6j > = (1 + n q )S qp , (55) 

with n q being a momentum distribution function. The averages <C 6^6^ ^> and <C ^> j n 
the case of normal effective vacuum are zero, while for a squeezed vacuum 

< b\b\ > = m 9 A(u; ? + u p - 2u s ) , (56) 

where the function m g is defined by the particular properties of a squeezed vacuum, and A (a;) 
is the discrete delta-function 



A(c) 



1 , u = 
, uj ^ . 
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The spectral function 7(0;) in Eq. (54) is assumed to be symmetric with respect to the central 
line, so that 

7(CJ S + UJ q ) = 7(CJ S - UJ q ) . 

Such an effective squeezed vacuum can be realized if atoms are inserted into a medium with 
squeezed collective excitations interacting with electromagnetic field. These excitations could 
be squeezed optical phonons or magnons. The interaction of the latter with photons results in 
the formation of squeezed polaritons. 

For the quantum field (54), the correlation function is 



«£ + W£(0» = £7 2 K) 



+ (l + n q )e 



-iU! q (t—t') 



+ 



(57) 



Calculating the quantum attenuation (49), we keep in mind that n q = n(u> q ) is real, while 
m q = m{uj q ) is, in general, complex, 

m{uj) = \m(oj)\e llfs . 

The final formulas will be simplified by remembering that \8\ = \uo — Q\ <C uo. Then substituting 
the correlation function (57) into Eq. (49), one finds 



n„ 



+ 



l + n 



with the last term being due to the squeezing, 

r s = -|m(ccO| 7 2 M 



— n{uj) 



7»r 

5 2 + T 2 



-17 



+ r s , 



(5? 



2 r cos ip s + 1uj s sin (p s _ Tt 

Auj 2 + r 2 e 



(59) 



The quantity (59) may be named the squeezing attenuation. 

Separating out of the quantum attenuation (58) its resonant part, one has 



n[uj) 



(60) 



As is seen, the resonant attenuation (60) is zero at t — 0. This means that, at the initial time, 
the main role in the quantum attenuation is played by its nonresonant parts, including the 
squeezing attenuation. 

To complete this section, it is useful to present explicit examples of n(u) and |m(u;)| in the 
case of a squeezed effective vacuum. The latter can be generated by a parametric oscillator 1 , 
with a squeezing field proportional to e cos(2cj s t + <p s ) . One employs the notation jj, = 7 — e and 
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v = 7 + £, where 7 is a cavity damping rate. Then, for a non-degenerate parametric oscillator, 
one has 



n{uj) 




(A s + k) 2 + ^ (A s -/t) 2 + /i 2 



1 1 



(A s + k) 2 + u 2 



1 



(A s - nf + z/ 2 



1 



m{uj)\ 




1 1 1 




1 
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(A s + K ) 2 + / u 2 (A s - K ) 2 + /i 2 (A s + /t) 2 + z/ : 



(A, - nf + z/ 2 



where A s = u — u s . The parameter k characterizes a two-mode squeezed field, representing the 
displacement from the central frequency of squeezing, where the two-mode squeezed vacuum is 
maximally squeezed. If the parametric oscillator is weakly nondegenerate, with \k\ <C |A s |, or 
degenerate, then 



By means of different parametric oscillators, it is possible to generate squeezed fields with 
various characteristics. 

7. Conclusion 

The method, presented in this report, makes it possible to combine the clarity of the semi- 
classical approximation with taking account of quantum effects. The latter are responsible for 
the appearance of a specific quantum attenuation. The vacuum electromagnetic field interact- 
ing with matter forms an effective quantum vacuum. Preparing different kinds of such vacua, 
one can regulate the value of the quantum attenuation. The latter, in its turn, can essentially 
influence the features of the following coherent processes. 

Notice that the squeezing attenuation (59) can change sign for the varying phase <p s . Thus, 
r s < for y? s = 0,7r/2. 

For preparing different types of effective vacua, one may incorporate resonant atoms into 
different media. There exists a rich variety of materials possessing various properties, which 
could be used for creating effective vacua with divers specific features. Just as a few cases 
of materials with interesting electromagnetic properties, one may mention biopolymers with 
photoexcited polar vibrations , photonic crystals with two-photon transitions , and plausible 
inverted superconductors 1 *^. 

The feasibility of influencing coherent processes by quantum effects opens a novel way of 
regulating the characteristics of collective atomic radiation. 
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